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ABSTRACT 


An analysis has been conducted of the capabilities of high 
pressure ratio, single-stage impulse turbines. Some recent experi- 
mental work has shown considerable promise for these machines. This 
work provides a tentative basis for better understanding the potential- 
ities and peculiarities of such high performance turbines. 

The present analysis indicates thaz well designed supersonic— 
for impulse turbines can be expeeted to produce static efficiencies 


closely approximating those predicted by well-known (subsonic) curves 


blade-speed 


of static efficiency versus ~—————————; - 
1Sencropr1e seus vecd 


Moparenivly.) there presemel exie comes tene neil sap lea ple 
method of blade row loss estimation in the case of supersonic flow. 
pmeextension to the Ainley-Mathieson loss prediction procedures to 
the supersonic case is proposed. 

Experimental and estimated supersonic-rotor loss coefficients 
Or the present report differ little in magnitude from well-established 
meonlts tor subsonic blades of similar deflections. 

TieSet eC. on EOrOr imeclaence. in Supersonic tlow, Nas 4 
profound influence on pressure and velocity distributions throughout 
the stage. The effect is inadequately explained by a recently proposed 
theory. The present analysis proposes an alternative method for esti- 
Mating the effect. 

A logical and straightforward procedure is developed for 


selection of the "optimum" design point for a single-stage impulse 


pyalleae 





turbine, given a required RPM, power output, and rotor blade root 
stress. 
The preliminary design of a large rocket-turbine is con- 
ducted as an exercise in application of the foregoing developments. 
This study was conducted as partial fulfillment of the 
requirements for the Degree of Aeronautical and Astronautical Engineer 


aoethe University of Michigan. 
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I. INTRODUCTION 


In recent years, considerable interest has centered around 
the development of high pressure ratio turbine stages for rocket- 
pump drive and other specialized applications. Desirable character- 
istics are reliability and light weight, with high efficiency as a 
@esarable objective of lesser importance. One type of machine attractive 
for these applications is the single-stage impulse turbine. This turbine 
combines low axial thrust on the rotor with potential for high pressure 
mab1os and specific work output. The high pressure ratios normally 
result in supersonic flow through the rotor, and very few experimental 
results concerning these stages are available. 

The present analysis is undertaken to explore the potential- 
meres of such stages with the benefit of some recent test results. the 
study is directed along three parallel, but distinct topics. These are: 
(1) analysis of flow patterns and loss distribution within the stage, 

(2) development of a scheme for selecting the design point to best 
achieve specified design goals, and (3) a typical design problem 


intended to illustrate the results of topics (1) and (2). 





TI. ANALYSIS OF LOSSES AND FLOW PATTERNS IN SUPERSONIC TURBINE STAGES 


1. General Discussion 

While a wide body of theoretical and experimental data are 
available (e.g. Reference 1) for predicting flow patterns and blade 
losses at low Mach numbers, little information is available at super- 
Semic velocities. As a consequence, turbines designed to operate in 
these regimes frequently fail to achieve the expected performance and 
fmessure distributions. Two examples are the turbines of Reference 2 
maemo woerein the flow patterns, pressure distributions, and efficiencies 
m@'rtered considerably from estimated design values. in both cases, the 
aeouloes would indicate efficiency was lower than necessary, had whe 
losses been accurately predicted in advance and compensated for in 
mae aesign geometry. 

The method of Reference 4+ was utilized by the N.A.S.A. to 
successfully correlate rotor blade momentum thickness and "specific 
blade loss" with "blade surface diffusion parameter" for one family 
Seeeoeveral transonic-rotor turbines. This method does not include 
Meemetrect of total rotor deflection angle, and has not been satis- 
itenerily generalized. The method predicts a grossly optimistic rotor 
fmementy coefficient for the turbine of Reference 2. The rover <detlec— 
tion angle is 140° in the latter rotor as compared to only 843° in the 
ieee .h. teamily, explaining the deviation. 

While the loss prediction methods of Reference 1 are compre- 


Memaive, they are only recommended by the authors in cases wherein 


22e 


-4- 


the flow is everywhere subsonic. Hence, an extension in generality 
to include supersonic rotor flow is essential to provide for accurate 
mew Channel layout and a more precise estimate of performance. 

Figure 1 presents the performance of the two supersonic 
meter turbines of Reference 2 and 3, together with that of the N.A.S.A. 
imensonic family of References 5 through 10. For comparison, there 
are the theoretical curves from References 11 and 12 and an unpublished 
curve from the Aerojet-General Corp. The theoretical curves presumably 
are intended to apply to subsonic stages only as no mention is made of 
Mach number levels. Slight deviations between the curves are explained 
by minor differences in the geometry and loss assumptions used in 
deriving them. They may be assumed to fairly approximate the range 
Of efficiencies to be expected for single stage impulse turbines 
Meamwines negligible Mach number effects. 

Comparison of these curves with the experimental results of 
mee high performance turbines shows agreement which is enlightening, 
if not surprising. It seems that no appreciable deterioration in 
Beeerormance may be attributed to Mach effects. The mediocre effic- 
moameics Of the turbines of References 2 and 3 are evidently due to very 
fen leaving losses, while the N.A.S.A. family has low leaving losses. 

ie erotic eGt ene ne rence. - slovouPehnernance tiated cite 
Bpuceuraging when compared to theoretical values. It is no surprise 
that this is the most recent and probably the most carefully designed 
of the three. This turbine is designed for high (140°) rotor deflec- 


en and maximum static efficiency. On the other hand, the N.A.S.A. 








Bamily of turbojet research turbines are designed for high total 
efficiencies, which average about 854%. 

Mes relative! yainterlor pertormancesen tne turbine sem 
Reference 4 appears to result from causes other than the high Mach 
mumiber level, although insufficient data are reported to provide for 
@ precise analysis of losses. The performance evidently is degraded 
byecxLremely severe flow separation from the rotor suction surface. 
This explanation is substantiated by the almost negligible change: 
mimtlow characteristics which resulted from a major modification of 
iene roLor profile. The modification consisted of reducing the blade 
thickness about 40% by removing metal from the suction surface. The 
negligible resulting flow change indicates the flow was separated in 
both cases, and demonstrates the ineffectiveness of profile modifi- 
Seulons downstream of the separation point. Since the profile 
feeteates tne criteria of Reference 14, as discussed in 1-5, consider- 
able separation is to be expected. In addition to the separation 
problem, the design apparently suffers somewhat from relatively 
large tip clearance and "lap". "Lap" refers to a sudden annulus 
eebireement immediately downstream of the nozzle exit. 

RevUE@haenL oO sens GeratLonseot cle sturaine wot Reference 2, 
Figure 1 also illustrates the margin of improvement - about two 
Pereentage points in static efficiency - obtained by increasing 
mevor flare enough to compensate for actual rotor losses. This 
mare increase allowed design impulse conditions to be closely 
eeroachned, and illustrates the effect of accurate loss estimates 


upon geometry and performance. Comparison of theoretical performance 





with the best performance of this turbine strongly indicates that 
any losses associated with supersonic effects must indeed be 
negligible. This deduction diametrically conflicts with standard 
works (e.g. Section 4.5 of Reference 1) which indicate a rapid 
increase in losses as the rotor flow begins to exceed local sonic 
velocity. Apparently a phenomenon exists which is similar to the 
"Sound barrier" drag increase of aircraft wings. Evidently, when 
the flow just exceeds local sonic velocity, strong and near-normal 
pmocks €Xist which cause separation. Once the flow is fully super= 
eeeee, the flow field appears to stabilize to the exvent that separa-— 
tion losses need be no more severe than in subsonic rotors of similar 
Beometry. The design criteria for minimizing rotor flow separation 
are considered in II-5. For the present, it seems reasonable to assume 
that the combined profile, secondary, and tip losses of the supersonic 
meLbor far outweigh any losses credited to supersonic effects. 
2. Development of a Method for Loss Prediction in Supersonic Turbine 
Blading 

inmeview vot Che foregoing examination of Mach mumber effects, 
the loss prediction methods of Reference 1] bear re-examination. ve 
basic Ainley-Mathieson loss parameter is 


y = (Total Inlet Pressure) - (Total Exit Pressure) 
(Total Exit Pressure) - (Static Exit Pressure) 


The total pressures are referred to relative and absolute velocities in 
the rotor and stator cases respectively. Following Reference 1, the 


woval joss coefficient for a blade row, Y ion GnecescUMrOt PilLer ILewl@cees: 


Je 





ue 


Yp» secondary losses, Y and tip losses, tee Profile losses for 


Son 
Seaver Or rotor blading of arbitrary reaction are empirical functions 
6f solidity, deflection angle, and thickness ratio. Profile losses 

are easily determined for blades of arbitrary reaction using the 
Methods of Reference 1 with Figures 2a and 2b, which are reproductions 
of Figures 4a and 4b of Reference 1. In the present case attention is 
mmeeeted Only to the special cases of nozzles and impulse rotor blades. 
In this case the value Yp is determined from Figures 2a or 2b, as 


eperopriate, for standard thickness + = .20O. The profile loss estimate 


merethe nozzle or impulse rotor blade in question is then simply 


= T/C aed 
"P= 9,50 ‘P (1) 


The secondary loss equation is developed from relations given 
in Reference 1. From Equation (4) of Reference 1, there is 


(Cp) sec S 
aS Remar rm a) 2, 
CT, 
From Appendix I, there is C; = 2(s/c) (tanB, - tanB, )cosB,, (3) 
and 


(C = gale) carey eee B, (4) 


eee sec 
where the vector mean angle, 


Je fat 


B, = tan [5 (tan B, - tan B,)] o 


m 


Here B is the gas angle-relative to moving rotor or stationary nozzle 


aS appropriate - and x and y refer to inlet and exit conditions 





respectively. Combining 2, 3, and 4, there results 


es 
per = mM cap (ean B,,-tan By) (6) 
i= scocriicient } iS an empirically derived coefficient which accounts 
femetne effect of acceleration on secondary losses. Figure 17 of 
Reference 1 correlates } as a function of hub-to-tip radius ratio 
and Ge Since these data were presumably gather for low sub- 


2 2 . , 
= (W,,/W,,) » where W is relative 


sonic velocities, it is assumed (A, /Ax) 
gas velocity. Based on the foregoing assumption, (i, Ww) should be an 
appropriate acceleration parameter independent of M. The assumed 
variation of \, with acceleration and hub-to-tip ratio is shown as 
meeure 4. It should be noted that experimental points of the original 
imeomre presented considerable scatter, whereas Figure 3 represents a 


mest fit of the data. 


The tip clearance loss equation is developed from 





(Cy) 2 
oy tae (7) 


of Reference 1, page 13. Equation (7) is attributed to Carter, and 
was developed on a theoretical basis for small deflections. However, 
Reference 1 points out satisfactory agreement between (7) and tip 
meets a2t large deflections, hence the extension seems justified 


in the absence of a better theory. Similar to (4) above, there is 


also 
2B 
(Je = a (Ey See (8) 
De Ke ee ; 


J 





SOs 


Upon combining (43), (7), and (8), there results 


oO 
(Ae 2 (=) cOsnaeen B,,)* (9) 


K cos ea 


For single stage stators, B. = O and (9) simplifies to 
x 


2 
(Male askil dale, 
Obes a 2(—) are (10) 


Similarly, tor impulse rotors with P,. = B, = B, there results 


4 


(Yx)R = 8(F)sin®p (11) 


The secondary loss relation (6) may also be simplified. 


In the case of impulse rotors (6) becomes 


COS NP me gage (12) 


imine case of nozzles, Figure 3 indicates that for any practical 


stator, 4. = -0055. Accordingly, (6) becomes 


Ss 
sinB si 2B 

(0G) aa) (OCIS) ie No Ce ees (13) 
Gospe cosB,, 


MieeStavor=cColLicicnbse calculated dn Ii-5.- tend to be 
excessively optimistic if (13) is used. Furthermore, the theoretical 
development based on Reference 1 and leading to (13) has so far made 
memallowance for effects of blade height on the magnitude of stator 
Beeconaary losses. However, the data of Figure 15b of Reference 1, as 
mem as Reterence 9 and earlier tests, confirms that stator secondary 
Hesses are inversely proportional to blade height, while solidity, 


meee eravio, and chord have negligible effect on stator losses. In 





order to resolve the apparent contradition, (13) has been modified 
feeecne coetficient 6/h, Cnosen Smpiricaliy based on test results 
mresented in II-%. Hence, the best estimate of stator secondary 
messes is felt to be 


2 
15 Sita 
—- —_ (14) 


oa, 
sec’s 
cos Ban 


The total losses of a blading row are determined by addition 


@eeune profile, secondary, and tip losses calculated above. That is 


Yn = Yo + Yoon + Xe aS) 


The Ainley-Mathieson loss parameter may be converted to a pressure 
@eetricient or velocity coefficient to simplify subsequent calcula- 
tions. 


One may write 


(ee 2 PP eT 
a ees eer a (16) 
2p jie p/P i 
YY y y 


(note P" becomes P' in the case of stationary blading). Then 


IL 


mieepressure coefficient is seen to be a function of the relative exit 
Mach number, Ms Soe volPeeas Clemo someOcrticto my Eante cas  pROperitcs. 


The conversion to a velocity coefficient is made using 


ae 


Cp = (1-7 M, a 





See 


which is developed in Appendix III of Reference 14. Solving for \y, 
there results 
y-1 


2 
We [2 ‘Gam? * - Cp 7 )] 


os 
= (18) 

Equations (17) and (18) were programmed on a 7090 computer 

and the results tabulated in Table I for y = 1.40 and ¥ = 1.25. ‘The 


results are also plotted for y = 1.40 and appear as Figures 4 and 5. 


4. Comparison of Loss Prediction Theory With Available Test Results 
The loss prediction theory of II-2 was compared with test 
results for 7 rotors and 5 stators of experimental turbines. These 
represent the currently available results that have been published 
with sufficient information to provide a comparison between predicted 
and experimental loss coefficients for blading rows. The comparisons 
covered a moderately wide spectrum of Y, My; derree Ol eteact lcm. 
deflection angle, solidity, aspect ratio and other parameters. While 
ithe volume of test data is quite limited, agreement between theory 
and experiment is reasonably good. The average deviation between 
predicted and experimentally determined coefficients is about one 
percentage point. The precision of both predicted and experimental 
coefficients is probably no better than plus or minus one point, for 
reasons discussed below. The comparative results are shown in Table II. 
It was necessary to apply certain assumptions and approxi- 
mations in compiling the data of Table II. ‘The more important of these, 


With reasons therefore, are: 





eae 


(a) The N.A.S.A. rotors of Reference 5 through 9 
were assumed to have tip clearance Ah/h = .02 
or Ah = .042", since clearance was not specified. 
The errors should be small, as Yx 5 Ym alg ste LiL 
cases. 

(b) Geometric and velocity characteristics are deter- 
mined for mean radius conditions, deemed repre- 
sentative of the stage. This assumes that three- 
dimensional effects are minor. 

(c) The effects of Reynolds Number variation are 
neglected. Reference 1 assumes (Re), =2x 107, 
while the experimental turbines operated at (Re), 
only 2 to 4% times this value. 

(d) The rotor pressure coefficients were not quoted 
explicitly in References 5 through 9. The co- 


efficients were deduced from 


n Pp, /P," 
(Cp)R = a 3 5 el es 
mm (CeCe) ng (Ga) ages 


which is developed in Appendix B. While this relation 
is exact, some interpolation was required to estimate 
Tm and Bee from the published turbine maps. The 
stator coefficients were deduced from Reference 9. 


4. Comparison of Typical Rotor Velocity Coefficient Estimates Obtained 
by Different Methods 


Figure 6 displays the estimated variation of rotor velocity 


coefficient, as a function of deflection, for typical impulse blading. 





Ae 


The subsonic reference curve is taken from Reference 12, and repre- 
sents a mean obtained from a large body of experimental results. 

The two supersonic curves were obtained following the method of 

II-2, for a typical geometry. The comparison illustrates several 
interesting features. The magnitude of the rotor coefficient estimate 
differs relatively little between the subsonic and supersonic curves. 
This is not too surprising in view of the comparisons shown in II-l. 
However, the shape of the supersonic curves at high deflections appears 
somewhat suspicious, as subsonic experience would anticipate a more 
rapid deterioration of velocity coefficient. Since there is insuffic- 
ient test data to confirm or deny the suspicions, one would probably 
be well advised to use the more conservative estimate, or perhaps 
pemiean. One also notes the increase of Vp with M, for fixed losses, 


Y. This trend is also seen in Figure 5. 


>. Detailed Analysis of the Effect of Rotor Profile Geometry on Losses 
Reference 14 reports extensive two-dimensional cascade tests 

of rotor profiles. These tests compared conventional constant- 

curvature steam turbine profiles to several experimental profiles 

designed to minimize flow separation in supersonic rotor channels. 

The inlet Mach number was 1.9 and deflection angle was 10°. 

Schlieren photographs and total pressure measurements indicated 

that the resulting "shockless" profiles essentially eliminated 

shock-induced separation from the suction surface. The over-all 


velocity coefficient correspondingly increased from .929 to .952. 





S132 


Reference 2 compared the conventional and shockless profiles 
in a test turbine, which also operated with several combinations of 


' and flare. Exact comparison between cascade and turbine 


rotor "lap' 
tests is difficult due to the limited number of configurations reported. 
Furthermore, the rotor inlet Mach number was about 1.6 instead of the 
design value of 1.9. The reduced Mach number resulted from off design 
nozzle pressure ratio, the reasons for which are discussed in Section 
II-6. 

the maximum rotor velocity coecriicient aehie vedm inne esc Goi. 
was .872. This peak performance was obtained using zero lap, but 
sufficient flare to approximately compensate for actual losses and 
approach the design exit flow conditions. In another test, the 
shockless profile increased total head efficiency about 2-1/2 percent- 
age points, under otherwise identical comlitions, while static efficiency 
changed little if any. These results will be discussed in more detail 
femesently. - 

The shortcomings of the conventional, constant-curvature 
passages are thoroughly explored on two-dimensional theoretical grounds 
in Reference 14. Essentially, it is demonstrated that constant-curva- 
ture passages cannot transform a uniform incident flow into a uniform 
curved flow. Rather, from considerations of a characteristics net, it 
is evident that the pressure on either surface undergoes a continuous 
zig-zag pattern of abrupt increases and decreases throughout the length 
Pe the curved channel. The pressure rises are expected to react with 
the thickening boundary layer and produce separation from the suction 


surface. Cascade tests demonstrated the predicted separation for the 


eonventional profiles. 





=e 


A second and complementary fault of the conventional profiles 
is the likelihood of excessive channel width caused by too large a 
spacing between blades. In supersonic rotor channels, the wider the 
channel, the larger becomes the uncompensated convex surface acceler- 
ation before the compression waves become incident. Hence, an 
excessively wide spacing produces an increasingly severe pattern 
of pressure fluctuations. Schlieren photographs also showed a 
definite tendency toward focusing of the compression "fan" if the 
Spacing were too wide. The abrupt pressure rise triggered separation 
end flow breakdown. It should be noted that supersonic rotor profiles 
are determined from channel flow considerations, and solidity is not 
an appropriate design variable. 

The "shockless" profiles of Reference 13 were designed to 
eliminate the foregoing objections. Their success in cascade tests 
was mentioned above. An example design is discussed and constructed 
in Section IV-5 and Figure 15. Essential features of the design scheme 
are: 

(a) Dual curvature, having an entry or transition 

curvature half as great as the central or mid- 

Chord curvauure. | [he transit lon eurnvauune 

extends to a position corresponding to intersections 
of the first characteristics originating from the 
opposite surface. Throughout the central region a 
uniform free vortex flow is theoretically established 
in planes parallel to the machine axis. A similar 


transition at the channel exit reconverts the flow 





abe 


GO uniiovmn parallel ext conGituons, 

(bo) A sufficiently small channel width to reduce to 
aeceptable Hevels the initial, uneompensated 
acceleration along the Convex suri ace.» “The 
maximum recommended width is that prescribed 


by a transition curvature having 


Yr - 
= > cos (sin7* = oem) (1) 
woul Gran es ae: 





(c) A channel width sufficiently small to preclude 
deceleration to sonic velocity along the concave 
surface. This requirement is generally less 
stringent than (b), unless the incident Mach 
number iS near unity. 

A comparison of cascade tests of the two profiles is seen 
in Figure {, where the data were collected from Figures 7 and 19 of 
Reference 14. Aside from the obvious general superiority of profile 
B, &@ closer analysis yields an insight into the nature of the compar- 
ative losses. The general spanwise variation of losses indicates 
the effects of secondary flows, wherein the end wall boundary layers 
are "centrifuged" toward mid-span by secondary vortices (see Art. 
14.7 of Reference 12 for complete discussion). A fair estimate of 
purely profile losses is probably obtained from the curve peaks, 
near root (and tip). The greatly increased losses in the central 
Span region point out the predominance of secondary losses -- in 


qualitative agreement with the estimates of Table II for the rotor 
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of Reference 2. The spanwise uniform superiority of profile B may 
be explained by a sharp decrease in profile losses with secondary 
losses essentially the same in both cases. Hence one concludes 
that the improvement shown by profile B is due to a reduction in 
profile losses caused by virtual elimination of shock-induced 
suction surface separation. 

The ratio of losses between the profiles, in cascade, is 
(Ya/Yp) cascade = 22 1.7, obtained from Figure 5 using the reported 
velocity coefficients. From Equation III-2-2, static efficiency is 
directly proportional to (1 + Vp) for impulse turbines with all other 


parameters assumed fixed. Thus, the relative improvement in yn due to 


profile B, if duplicated in the turbine, would be a. factor of only 


1+ 922 | One: 
1 + .929 


The reported data of Reference 2 allows at least a rough 
Comparison of losses obtained using profiles A and B in the test 
turbine. Total efficiencies were .610 and .635 respectively, with 
only the rotor profiles changed. From Equation II-3-1, (Cp)p/ (Cp), = 
1.06 for stator velocity coefficient and total pressure ratio fixed. 
Also, using Figure 4 for My ~ 1.4, together with an estimated Y, = .84 


from Table II, one may show (Y,/Y¥ ~ 1.18. Using Figure 5, the 


3 /Rotor 
improvement from, say Y = .64 to Y = .84/1.18 = .71 is equivalent to 
about .013 increase in Vp. Because n is proportional to (1 + VR)» 

it is understandable that the static efficiency improvement, due to 
the improved rotor profile, was indeed negligible in the test results. 


several important observations may be made from the results 


of the cascade and turbine tests of the two profiles. It is obvious 
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that the design features of the shockless profile B prove beneficial 
in all cases. The relative improvement of profile B is donsiderably 
less in the turbine, however, under this single set of test conditions. 
This result, together with a velocity coefficient reduction from 
Voascade = ° 922 tO Vp = .872, may be at least partially explained by: 

(a) the likelihood of moderate flow separation in the 
turbine rotor, while vircuaily sene  oceurred. in 
cascade tests of profile B 

(bo) the degrading effect of tip losses, not present 
in cascade, which also are proportionally more 
detrimental to rotor profile ss. 

(c) off-design operation of the rotor, which 
operated at My, = 1.6 instead of 1.9. 

(a) other three-dimensional effects, such as radial 
pressure gradients, which would be expected to 
degrade the performance of moving blades. 

Tn addiGion te thie elleccus sol sprouse —ecudnpe ened erence: 
also illustrates another important facet of rotor losses and the 
Mecessity for their accurate estimation during design. In one test 
series, an improvement from 7 = .50 to .52 and Wp = .844 to .872 was 
accomplished by a slight increase in the rotor flare. Both tests 
used the "shockless" profile. In the first case the flare was based 
on an estimated Vp = .91 and, due to the magnitude of the actual 
losses, the exit swirl angle B), was constrained by continuity to be 
considerably less than its design value. The increased flare approx- 


imated actual losses and produced the increased exit swirl and efficiency. 
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items evident that an accurate estimacre of the actual velocity 
coefficient may be considerably more important than the question- 
apie ability to increase it a few percentage points. 
6. Analysis of the Effects of Rotor Incidence on the Supersonic Flow 
Through Turbine Blading 

On the basis of experiment, Reference 14 concludes that rotor 
blading with non-zero edge thickness, in supersonic flow having a sub- 
sonic axial component, is constrained by continuity to operate at a 
Unigue incidence relative to the oncoming gas flow. This incidence 
@ezermines the gas direction relative to the rotor. Likewise, the 
absolute discharge angle of the stator, as well as effective stator 
meieearea, is determined uniquely and is virtually independent of the 
mater geometry. The rotor incidence effect, if unceompensated, can 
drastically reduce the stator discharge area, pressure ratio, and 
Mach number. The stator pressure ratio of the test turbine of 
ieterence 2 was reduced to between 1/2 and 2/3 of desiga value, for 
Beep le, due to existence of rotor incidence of hee: ital seyident 
fore must be taken to estimate the rotor incidence and thus determine 
mieerequired blading geometry to minimize losses and obtain the 
@eaired pressure and velocity distributions. 

References 2 and 14 set forth a theory for determining the 
Unique incidence, i = Bo - Ba. This theory will be stated briefly, 
Smeets validity examined, in view of the importance of the incidence. 

The incidence is postulated to be that for which the expan- 
sion, corresponding to a Prandtl-Meyer turn of Bo - Bz, is exactly 


equal to the geometrically available expansion ratio (s cos Bz-t )/scos B,. 
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There is assumed to be a slender wedge leading edge with attached 
shock, and negligible total pressure loss accompanying the expansion. 
fenee, the problem as postulated resolves into simultaneous solution 
er the continuity and Prandtl-Meyer relations. 

The Prandtl-Meyer relation in this case simply v + B = 


Constant or 


v(Mys) = Vii.) + (Bp -B5) = vy.) + 4 (1) 





A convenient form of the one-dimensional, steady continuity equation 


may be developed form 


is 
0 VA = (S7) (Mt VyYRT) A = constant (2) 
hence , 
nes constant (4) 


But T = se and P= Bs 
iL dk Je ma G i = M27) Y-1 


Since T" is constant through a constant diameter rotor, 


PU MyA 1) 
= = Weenetant ( 
= M,?) /2( 
Examining stations (see Figure 9) 2 and 4, just upstream 
and downstream of the leading edge, one writes 
pat 
pe 


-1| . 
ie a+ Mas? Mas 7/ ae) (5) 
> 1" -| 
Ay Eo Me ees 7 Mug? 
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From the blading geometry, there is also 


A cos B. - t/s 
————— (6) 
Ay COS Bo 


Using Cy = Pz"/P5" and combining (5) and (6), there results 


y-L 
= cos Bz - cea? "Wo 1 + Mas y+1/2(y-1) (7) 
Cp Mar ee eM, oO 
2 


Equations (1) and (7) may be solved for two unknowns. 
Typically, one knows M5» Ba, and t/s and wishes to find Mir and Bo. 
The incidence is thus determined for a specified incident Mach 
Mumber and rotor blade proiile. 

Reference 14 illustrates the theory for an example rotor 
having Bz = 70° and May =1.9. Reference 2 compared estimated and 
experimentally deduced incidence, with fair agreement, for Bz = 70° 
and Mi. = 1.61. See Figure 8. Isentropic conditions are assumed 
Throughout. 

Upen close Seruviny, elOvever. wie Oo veve Uncen may. be 
shown to become absurd as the rotor blade angle, Bz, doe Geuuced 
somewhat. The anomaly is easily demonstrated by a simple numerical 
example. 

Suppose the (relative) incident flow occurs at Mir = 166: 
The corresponding Prandtl-Meyer angle for y = 1.4 is Vo = eon. 
Further suppose a fixed incidence of any reasonable magnitude, say 
3°. Now using Equations (1) and (7), calculate the allowable edge 


thickness ey he for various Bz. Assume isentropic conditions, 


Cy = 1.0. 
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By (1), vz(M,_) = 17.9° which yields a corresponding 
oS 
Mur = 1.703. Equation (7) becomes 
60 ie peOle one |e 
t/s = cos Bz - cos (Bz + 3°) oe oc i ea 
? ? 1.705 |1 4+ .2(1.60)4 


= cos Bz, = 1.0/5" cos (Bz + a) (8) 


For Bz = “(Okan t/s = .02/. However, suponmeuoauce tne Ba it is seen that 
t/s decreases monotonically until t/s = OF for Bs - 52°. The same 
conclusions are reached regardless of the selected numerical values. 
In particular, it is evident that breakdown occurs at even larger 


eat Gy = 1-0 or My pete 


2) 

The physical explanation of the breakdown of the theory is 
that, at sufficiently small blade angles, the turned and accelerated 
Prandtl-Meyer flow (assumed) fails to encounter the requisite geometric 
expansion necessary to satisfy continuity. In addition to the fore- 
Soing arguments, it is also evident that the isentropic, irrotational 
Prandtl-Meyer flow assumptions cannot hold in an actual turbine rotor. 
In article 9.10 of Reference 12, Vavra has pointed out the inherent 
danger in applying two dimensional analysis to such three dimensional 
flows. 

Aside trom the merits of the foregoing theory, rotor incidence 

remains a physical reality with very important implications. Hence, a 
more satisfactory explanation of its nature must be sought out. 
One must agree the continuity relation (7) is valid. Having 


Gecided the Prandtl-Meyer relation is not applicable in this case, 


there remains one more unknown than available equations. 
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Expanding the Mach number term of (7) in a binomial 


series, it may be shown that 


Y-l  -2\ y+1/2(y-1) 
Mp Le Mis =lite  M, - 1) + O(e) +. 
M 7-1 
me ae a Mio 
where ¢€ = “rs a De 


Substituting in (7), there results 


Ytl ye 
a a rene Bo (9) 


Cp 


t/S = cos Ps a cae e( 


From (9) it is evident that for moderate incident Mach 
numbers, the magnitude of the incidence is relatively insensitive to 
emell! differences between Mur and Mio: 

Equation (9) was compared with the experimentally deduced 
Mmeerdence of Reference 2. For design speed conditions, one is abile 


ice depuce ch. + 1.9) = Oh. Assuming ¢€ > 0, one then concludes that 
oh <Cp<1.00 ; O'< ele ome 


As stated earlier, a second valid equation in the unknowns 
is apparently lacking. Since it seems probable € << 1 in cases of 
interest, one might assume ¢€ = 0 and rewrite Equation (9) as 


cos Bo 
Cp 


t 
See CS ee (10) 


The theoretical incidence of Equation (10) and of Reference 
14 was compared with the experimental incidence of Reference 2. The 


comparison is shown in Figure 8. There is fair agreement, at the 
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design point, between the theoretical estimates for isentropic 
conditions. Furthermore, Equation (10) is seen to agree with 
experiment for Cp = .o4 at design speed and .90 at zero speed. 

This trend seems reasonable and agrees with the expectations of 

the authors of Reference 2. Higher relative losses at design speed 
are attributed to the reduced Mach number with correspondingly 
stronger, and possibly detached, shocks. 

By examining the design speed case further, one may deduce 
mic apparent proportion of total rotor loss, ocecurrings at the leading 
edge, which results in the 4.5° observed incidence. Over-all rotor 
losses were reported as Vp = .872. Using Figure 5, this loss corre- 
sponds to (Yp)p = 1.12. Similarly, from Figure 4 the apparent lead- 
ing edge loss, Cp = .84, corresponds to (Yrn)R = .25. Hence, the 
apparent ratio of leading edge to total rotor ioss is Beiey ei. = 2/9. 
As noted above, this apparent ratio represents an upper limit of 
actual leading edge losses and, in effect, incorporates the effect 
Bimmany flow acceleration, ¢€, which may occur at the rotor entrance. 
From Equation (10), for a given blade (or gas) angle and incidence, 
leading edge losses would be expected to increase with t/s. The 
foregoing considerations are utilized to estimate incidence in the 
design problem of IV-5. 

ihe result of a single experimentemardhy conssi vubes 4 
reliable design rule. However, there is provided a reasonably con- 
Sistent account of the interaction between gas angles, blade angles, 


and rotor geometry. As mentioned previously, care must be taken to 
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aecurately estimate incidence and preclude serious rotor-stator 
mismatching with resulting losses. One test cited in Reference 2 
clearly indicates the criticality of the problem. Rotor blades 
having negligible edge thickness were run. During the course of 
the test, the nozzle pressure ratio decreased from 16/1 to 13/1, 
corresponding to the accumulation of only a thin carbon deposit 
on the rotor blades. Nevertheless, the buildup caused the 
incidence to evidently increase some 2° in order to produce the 


observed effect. 





III. A RATIONAL TECHNIQUE FOR DESIGN POINT 
SELECTION FOR SINGLE-STAGE 
IMPULSE TURBINES 

ieee General Discussion 

Given a design goal in Germs of skit) and power required, 
expansion ratio, working fluid, and thermodynamic inlet properties, 
one is faced with the task of selecting a design point to satisfy 
micerequirements. Literally an infinite number of design points, in 
werms Of Dilade speed and gas angles, can meet the requirementes. hence, 
the designer must choose the most desirable combination of high 
efficiency, compact size and weight, and tolerable stress limits. 
An obvious approach is a trial and error process. This results in 
laborious calculations of debatable final success in determining 
the "best" design point. 

A single-stage turbine, especially one with a high expansion 
Tatio, will tend to have sizeable leaving losses. Hence, efficiency 
mopoetermined primarily by blade to isentropic jet ratio as seen in 
Figure 1, while the blade speed is limited by material stress con- 
Siderations. Experience and the mission requirements generally 
dictate desirable blade materials and allowable stresses. Thus one 
teat pe well advised to restrict the infinite number of possible 
design points to only those which satisfy some specified blade root 
memess. One could then proceed to decide upon a desirable size, 
weight, and efficiency combination to meet design requirements. 

The following method is one way of logically deciding upon 


the "best" design point. It requires neither machine calculations 


Ble 





nor infinite human patience. The technique is applicable to single- 
stage impulse turbines, totally independent of Mach number levels. 
The loss coefficients which enter may be estimated by any means what- 
soever, and are not related necessarily to loss prediction methods of 
II-2. The nomenclature of the development is defined in Appendix A. 
The method is applied to a design problem in IV. 
Beeovacic Efficiency as a Function of Velocity Coefficients and Rotor 

et lection 

In terms of Figure 9, Euler's Turbine Equation may be 

expressed as 


Co 


ber symmetrical impulse blading of the present case, 
Bee alee es eee 


Equation (1) may be rewritten, 


7 = 2 oT Gl Vp) Wo sin B, 
O 
or 
U Wo 
ee ee ee Vp) sin B (2) 
Ce 


By definition, ¥_ = V,/C_, and (2) may be expressed 
s rox Se 


2 
Ss 


(1+) 06 sin B (3) 
From the geometry of Figure 9, 


Q@ = Vl -£ coséB - Ssin B (1) 
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Hence (4) may be expressed 
n = ae (1 +¥,) sin B &W1 - 6“cos°B - SsinB) (5) 


Vp is a known function of both B and § , in general. ¥. is inde- 
pendent of §, but weakly dependent on B - through @. Equation (5) 
Mere ceecalenlation of efficiency for selected design point coordinates 

& , B. While Equation (5) might be considered to "organize" the 
search for a Suitable design point, it places no restriction on the 
blade stresses. 

Incidentally, Equation (5) also illustrates the relative 

mmportance of Ve as compared to Ve For example, a one percenuace 
Peint improvement in Ve should be reuphly Equivalent te aaneur pornic 


improvement in VR> since Va ~ Vp ® ds 


pee Limiting Stress Requirement 

For a selécted blading material of | specified) taper tauie, 
including constant cross section blading as a special case, one may 
emiew that blade root stress is proportional to the product of radial 
meeeleration and blade height. if bending stresses are neglected, 


2 
Go ~ hh a (6) 


by CONnCGinuity, the mass flow rave is, 
O 
fee oecnsity x axial velocity x annular area, or 
O 
W~ Q),W), cos 2 hi me 


the exit density, depends on pressure and temperature at 


Py 


mee blade exit. The pressure may be assumed independent of the vari- 


ables and equal to exhaust pressure, provided the axial component of 





ele 


velocity is subsonic. The temperature will vary somewhat with 


n=7 (6,8). If one neglects this variation, Equation (7) becomes 


O 
W ~ W), cos B Rh (8 ) 
; dk 
For fixed power output, W ~ z* Since Wy = VpWos Equation 
(8) becomes 
~ Wp Wo cos B Rh (9) 
Substituting, 
a eee 
NVR WoR cos B 
into (6), there results 
U\2 + 
one (=) (10) 


R Wr Wo Cos B 


One notes that U/R; the angular velocity, is fixed in the 
Beesent case. Further, letting Wee S Vz ~ § , by assuming V3 


independent of § and&, Equation (10) becomes 


Se et 
S NWR cos B 


or 


a) 


lll 
to 
: 
: 


on’R 8 cos B 


The parameter B has stress units and may be viewed as an invariant 
property of the particular design specifications and working fluid. 


Experience indicates that B is indeed virtually constant over quite 
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a sizeable range of §, B variation, as illustrated in IV. 


If one now selects a specific working stress 0 as a design 


In 


objective, there results D = = Ne ® cos B (1) 

o 

mis, tor a specified stress, the parameter DP is 4 unique constant. 
Furthermore, any design point satisfying (12) automatically satisfies 


the imposed stress requirement, ¢G. 


4. Combination of Performance and Stress Requirements 
By combining the results of Equation (5) and (12) , 


2 : iqieenenecy : = D 
2y~(1 + Vp) sinB $(N1 - cos<B - $sin B) * 7 bcos B (13) 


or 


: D 
Ger i -§© cos®B - §sin B) eV. Wala + Vp) sin B cos B- (14) 





For convenience, one may define parameters 


D 
Perce Dl - 
aN VRQ + Vp) sin B cos B (ha) 
rug Seen - $-cos“B - §sin B) (14) 


Note that F" is a universal function of 8 and B and could be tabulated. 
Since Wp = Vp §,6) and Y, is assumed independent of B, while 

D is a known constant (for a specified 7), one may evaluate Equation 

(14). There remains an unlimited number of admissible § , B coordinates 

which are Known to satisfy the imposed stress limit. The problem is 

the selection of the most attractive 6, B that represents a final 


@esign point decision. 


As an aid in this decision one might examine some variable 
additional relationships. 


For example, from (13), 


7 D 
n(§,B) = t= 6coERe ely 


The rotor relative Mach number variation may be examined by noting, 








We J RT V7RT5 


Meee mermore, if Vo = VC, is independent of B, as previously assumed, 
it follows that In is likewise constant. Hence one notes, 
Mo~ &§ (16) 
2 
Variation in radius is simply determined from 
U ~~ RR» @o a 
Reference 15 gives a useful empirical estimate for single-stage rocket- 


turbine weight as 


m,lb. = 7On (R,ft.)© = 1.53 (R,in.)? (18) 


finally, since o ~ Kh ~ Oh, for constant RPM, blade height variation 


Mey pe deduced from 
—— = constant (19) 
De Frocedure for Application 


Given power and RPM requirements, together with available 


expansion ratio and turbine inlet conditions for a specified working 
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fluid, one may proceed in an orderly manner to select the "best" 


design point, as follows: 


(A) 


(B) 


(C) 
(D) 


(E) 


Based on experience and operational objectives, 
one makes a tentative selection of blade material 
and one or more values of allowable centrifugal 
blade root stress. 

A trial design point is selected in terms of an 
arbitrary velocity triangle. As a minimum, one 


evaluates, for the trial design point: 


7 VR 
O 

W B 

oO Mr 
R Mary 
h o 
¥, 868 


A single calculation of this type should suffice 
provided the calculated stress is within, say + 50 
percent of the selected values, Gc. 

Calculate the invariant B - Equation (III-4-11). 
Calculate D corresponding to selected o - Equation 
(III-4-12). 

Plot F' against B (for a reasonable regime of 8B) - 
Equation III-4-l4a. Consider Y, constant as calcu- 
lated in step (b). Consider Ne = Vy (BM, ) ~ Vp (B); 


using constant Mary, as calculated in step (b). The 





(F) 


(G) 


(H) 


(I) 


a2 


curves F' thenrepresent constant stress lines. 
Slelstealgabeabhg=| 10 cheno (cial ieardsihediie - 

Plot F" against B, for selected ® = constant - Equation 
TI-4-14b. Intersections of the constant § and constant 
o curves determine unlimited compatible (§, B, oc) 
design point possibilities. 


For any desired (&, B, co), calculate: 


7 - Equation (111-"-15) 
Mar, - Equation (III-4-16) 
9 - Equation (III-4-4) 

R - Equation (III-4-17) 
U - Equation (III-4-17) 
m - Equation (III-4-18) 
h - Equation (III-4-19) 


As an additional aid in visualizing the variations, 

plot n against B for o = constant. The addition of 

$= constant and 9 = constant points provides for ready 
visualization of Ehe variation. ine, M5» QO, R, U, m, and 
h along the constant stress curve. See Figure 11. 

From the plot of (H), one may then select the most 
attractive Combination et  earures te fix tne mesien 


point (&, B, o). Estimate n, M,_, 0, R, U, m, and h 


3? 
for the selected design point. Repeat (B) to verify 

the estimates above. One notes that the approximations 
introduced into the technique are now discarded. A 

second iteration of (C) through (I) would appear unnecessary 


unless (1) indicates intolerable lack of agreement between 


the rough and smooth design point estimates. 





IV. ILLUSTRATIVE DESIGN PROBLEM 


i specifications and General Design Philosophy 

An example of a preliminary design is presented to illustrate 
the application of previously discussed methods of loss prediction 
and design point selection. The specifications and gas properties 
arbitrarily chosen are those which might correspond to turbodrive 
mequirements, of a large liquid propellant rocket using reasonably 
energetic chemical propellants. 

Assumed output requirements are: 28,800 SHP at 15,000 RPM. 


Assumed gas properties are: 


I 


molecular weight, ™m 129 1b/mole 


ieee 


Speci1 1c Meay fano.s 


constant pressure 


specific heat oa soe] Btu/1b-"R 
total inlet pressure P!' = 450 psia 
O 
exhaust static pres- 
sure Py = 40 psia 
ze O 
gas constant, R = 128.8 ft-1b/1lb- R 
total inlet temp- T! = B2200°R 
eracure, 2 


The turbine is to be of the single-stage, impulse type to 
minimize axial thrust on the rotor bearings. A reasonable compromise 
between static efficiency, size and weight, and simplicity is to be 


achieved. 


~44~- 
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2. Design Point Selection 

The design point selection method described in Section III 
aseco be followed. A single, arbitrary trial design poin¥tserves 
as a basis for selection of the final design point. 


Tae available isentropic enerea, els, 


Py y-1/¥ 30, .20 Btu 
—= ! = omens = ‘ — —_—ea 7 — re 
Altg = Cp T lal ea ] 827 x 2200 [1 To) ] = 760 a 
O 
The isentropic jet velocity, 
C= V 2g JAH. = ¥5 x10" x 760 = 6160 ft/sec 


The isentropic Mach number, 


(M Payee! 3) Ne Gee ye te 


) = 
2/IS aS 
yl Py 


The subsequent steps correspond to those of Section III-5. 


STEP A 
Suppose consideration be given to blade root stresses of 
30,000, 37,500, and 45,000 psia. Most suitable blading materials 


are known to have densities of about 0.30 Lb/in?. 


euuP B 
For the first trial arbitrarily select a, = 65°, U/C, = .3. 
From Figure 1, the expected efficiency is .60 < 4 < .67. The mean 


blade speed is, U = .3 x 6160 = 1848 ft/sec. The mean blade radius is 


_ 2846 x60) , io = ihe" 


2x xX 15000 





Zaps 


(1) Pelemiuneuion Of —olauor Velocity Ceemicient.. Vy .— 


SE et ce a ap ere ire eine Biel Abe a oi ae ee ee ES 


The mass flow rate is 
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* Soe re. eo cao 
The annulus area, 
O O 
A = ernRh = ene — W 
2 2 P 
PoVax 2 ey Car eosd: 
z RT SS) a 
oe 
WRT» 
h = 
Vg Po Co cosy 2nR 
For an estimated , x 95. 
AH 
is = ts 21S . 2000 - .95% 16 - 1370°R 
_ = Eee 827 
and 
1b EE ub O 
: Nera leo Tpoop x 1370 R oom 
eae” 195 x 30 ae 6160 2. eos 65° x On x 14.1 in === 3 
in sec 


approximate annulus height. 


The nozzle secondary losses are calculated from Equation (II-2-14), 


ales. sin“ 
Gage le a 
2 cos ome 
where 
2 ten tana.) = tant (22) - 47° 
ont 2 2 a 2 
Thus 
» ef cO 
cae = lS esinsoo - ah 


sec’S 1.12 cos 47° a 
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Stator profile losses are estimated from Figure 2a to be Cae = .06. 
Total stator losses, (Yo), = (Yoecle + Cole < .20. From Table I-d, 


for Mp = 2.2, the velocity coefficient estimate is y, = O74. 


(2) Determination of some velocity triangle components - See Figure 9. 


Nig SO -974 x 6160 = 6000 ft/sec 





VeseD = Vy, cos A, = 6000 cos 65° = 2535 
- Vo sin A - U - 
By = tan"+ (=——=$——) 5 Bet Ee Blue 
aX,e2 25 35 
Wo = —82 - 2535 _ = 4390 
OT oie cos 54.7° 
$ = Wo/Vo = 4390/6000 = .732 
9 = U/Vo = 1848/6000 = .308 
in ae! - y,° SIS = 2200 - .97K* 760 _ 1 506°R 
Cp 827 aaa 
. 
Revised M. = a= Esudo. = 2.29 
N yeRT, 72.0 N1328 
ee = 4390 = 1.675 


72.0 Ts 2620 


(3) Determination of rotor velocity coefficient, Wp 


Rotor secondary losses are, using Equation (II-2-12), and 


estimating \ = .03 from Figure 3, 


CRar = 48 gee = 48 sino pice a = ae 


Tip losses, assuming Ah/h = .03, from Equation (II-2-11) 





= 


ae = 8.0 = sin'B, =o Rone OS sin® 54. 7° = 5 ike 


Profile losses are estimated from Figure 2b, (Yp)p x .12 
Estimated total rotor losses are (Yn)p = .60. 
0) Bh Sensbienes| welll Gndve Wen ADEN Ne 1.6, the estimated 


velocity coefficient from Table I-d is Wp = .90h. 








The rotor coefficient compares with an estimated Wp = 898 


from the subsonic reference curve of Figure 6 for B = 54.7° . 


(4) Completion of velocity triangle and performance calculations 


W, = VpWo = .904 x 4390 = 3970 





B, = Bo = 54.7 





ae = Wy, cos hy, = 2290 
Tes Nee 0 + (Wy, cos By See Neumeoce 
, ax, 4 ed alt paths 
V 
-l tax 4y _ -1 (2290) _ O 
A = cos V? cos E80" Mes 


. ; . 
work output, W = cal - We (1 +V¥R)sinBo _ L848 x 4390x 1.904 sin 54.7 
gd BJ 2 ee 


W = 504 Btu/lb. 


2 20 ee 
(E-aiCO), Maen 


peeolc efficiency, jy = 
This estimate appears reasonable compared with Figure 1l. 


(5) Temperatures, annulus (rotor blade) heights, and rotor blade stress 
Since total relative temperature is constant through the rotor 


of constant mean radius, 








rT," = T Woo = fe ee 
© 2gJCp 2gICp 
S2nee 
ae Won 2 3. i390? 2 
ies aie? k= *p 2gJC GS Vie) = 8 "LS & 16 ee 
fp) pea culecgeaa Mm = th =) a 
: 4 720 V@, a 
2 
V Yn 
ae. ae 2680 : O 
et), 6 6= «Ty C+ 2EICp HLS ert. = 1589 R 
Check work output, W = cp AT' = .827 (2200 - 1589) = 505 Btu/1b. 


fits Constitutes a check for consistency of the velocity and temperature 


fascribuLions. 


The mass flow, 


220 


O 77g SHP 
el OS eS SIO Th lb/sec 
W 504 Sa 
O —_ 
: Wee 
the annulus, or rotor leading edger eign ic Hy = 
enRPoVax 2 
NG) x Hele <6. x 1328 ‘ 
ae oS ee Sees 
ele re INL So SO ne Gs — 
WR 1413 ab 
Ay ~ pet) =hy tel = EOE os aa a ee 
2n RP) V Vy 32 : See 
oS ST weal To YR » 


The blade root stress is calculated assuming a blade of height 
hy and density . 30 lb/in?. The centrifugal stresses only are considered, 
as bending stresses are expected to represent a small fraction of total 


root stress. 
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h Ea 
The root force, F = ma = (pA, a (RL oF ft/sec“) 


hy R_ (2x 15000)* 


r= ie 
g 12 60 
hy, R 30000 xe : 
co = P/A., = i eh ee) = 1910 h)R psi 
ee 12) 60 ee 
for h, = 1.207" and R = TRS: =) Gof elcial 
(6) Summary of step (B) results 
7 = -663 y = .904 


oe Mel Walace 9 2 _ oh 


o = 32600 psi My = 1.675 
RV a May Se hs 5) 
i a ak elo Grr Wied 
Ve = -974 Se 2 
SleP C 
B = ONVp & cos B = 8250 psi , the invariant 
Slee D 
(1) B 8250 a5 
D = = =) ae 
=) 30000 ace 
2 cA 
™ es . 220 EES, 
OTEP E 
Constant stress lines represented by F! \+ = eRe a 
2V5VpR(1t+ Vp)sinBeosB 
are plotted in Figure 10. For this calculation, y. = .974 = Constant, and 


Ss 


Ve = ¥R(B, Man, = Vp(B, 1.45) is assumed. The total rotor losses are 





Sie 


esimated following step (B-3) to be 


(Y (Y + ¥ + (Y 


m2 
ae ele (Cn eneo tele 


DR = p)R i 


The velocity coefficient, Vp> is )then casily secleileted se ase ins ote 
(B-3) as a function of variable B with Mon = 1.45. A simple working 
mos of Vp = Vp (B) constructed from Table I-d is a convenience in 


Steps (E) through (H). 


SleP iF 
Curves of F"(f) for i constant are plotted in Figure 10. 
The range .52 < 6 < .76 represents practical limits of interest and 


meerements of .O2 are satisfactory. 


STEPS G AND H 
From intersections of $= constant and o = constant curves of 
Figure 10, Equation (III-4-15) was solved for yn. Plots of 7(B)_ 

6 =const. 
are shown for the selected stresses in Figure 11. While not expected 
to be quantitatively precise, Figure 11 condenses the -design variables 
into an easily digested presentation. Since n, Mir? ley UE hy, depend 
wueyeon 9, & , B, the selection of a deSign point to fulfill 4 


Specified mission is considerably simplified. 


im che present case, one Caniwirite. 


Bee U-416), a 1.675) eee 
(Eq ) oo (Pas 29 § 


(Eq. III-4-4), 9 = 9(6§,B), plotted for convenience in the 


WSVelMNS (Ene Wo wees sh. 


o 
(Eq. III-4-17), R w 14.1 es = 45.8 0 
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epeLll-l-17 ), is) ele 6 ee eGonone 
. 308 ae 
(Eq. III-4-18), Weel (ls ee = 3220 35, ile 
(Eq. III-4-19), hy ® ages % x 16207. = 42 lene o 
9 32600 2 eee 


STEP I 
A working stress must be selected. The blade temperatures 


may be estimated, for a (turbulent) recovery factor of about .9, to 


De, 
a ee ee 
blade - “2 °°” 2gJCp ? ha3xi0t) = taal 


Prom Figure 11, one notes that efficiencies are reasonably geod tor 

o = 37500, with relatively little gain to be achieved at higher stresses. 
Experience would indicate o = 37500 to be a logical choice for a 

rocket turbine application of short running time. Ample safety factor 
should be available for a blade material such as Inconel 703C at this 
imemoecCracure. 

In Figure 11 for ¢ = 37500, the regime, .60 < § < .68 appears 
most practical. A rough idea of the trade-off between fuel consumption 
and turbine weight may be quickly made if desired. For example, com- 
paring points at $, = .66 and So = .62 shows the latter turbine would 
weigh approximately An = 3220(0@ - 08) = 3220(.4a" - .37°) = 127 lb. more. 


The corresponding reduction in fuel consumption would be 


0 O 550) GieG Oba 1 .707 x 28800 a i 
AW = Wo - Wi = —— =) = le oa 


1.34 1b/sec. 





elas 


The choice of O= .62 would appear advantageous only. if the burning 
iEnone ceeds 


om . tet =. 95 see 


A 1.34 


provided the objective is simply minimizing the sum of turbine weight 
and total fuel consumed. 
In the light of foregoing considerations a choice was made, 


selecting: 6= .64, 68 = 60° with the assurance o » 37500. From 


Figure 11 and Step H, the additional estimates are immediately available: 


QO ~ .395 

i eS 

O 

W * 34 1b/sec 
Ma ~ iL lNe 

ee koe 

aa eo 

U * 2370 ft/sec 


DOO a 


= 
22 


The design point selection technique has allowed a logical 
and systematic choice of operating design point. It appears that a 
considerable improvement over the original arbitrary choice has been 
achieved. However, detailed design calculations must follow to verify 


the choice and remove the inherent approximations. 


4, Detailed Design Point Calculations 
The absolute gas angle required is 


O4 


> = cos7+ ( § cos B,) = cos7+(.64 cos 60°) = volo 
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and 9 =1.0 sin M - §cos Bo = sin 71.35° - .64 sin 60° = .392 


(a) Recalculation of stator velocity coefficient, ¥., 


In the manner of Step (B-1), ¥, must be changed to incorp- 


orate the revised M and nozzle height. For an estimated ¥, = .%, 


To = 2200 - 96° _760 = 1355°R 
627 
pile been re Mop) ‘i ; 
h See, a ximate 
« 20 (.96 x 6160 cos 71.35°) On x 18.1 ee 
ennulus height. For On, = tan? (5 tan 71. 35°) = 55) , 


iB eatin Jal, So" 


(yo es A 
(920 COss oo 7 


ne) > 22k (Yb), =~ .06 


Total stator losses, (Y,,4), = -28. For M, ~ 2.2, Table I-d gives 
Vg = +965. 

It is evident that both (Cp), and ¥, depend on the, as yet 
unknown, value of M, - estimated to be 2.2. In subsequent calculations 
it is assumed ¥, = .965 to determine the actual My. The value of Cp 


will then be revised, where necessary, to coincide with the newly 


Calculated Mp. 


(b) Determination of some velocity triangle components 


Vp = .965 x 6160 = 5950 ft/sec Vax,2 = 2950 cos 71.35. = 1900 
Wo = 8Vo = .64 x 5950 = 3810 U =W, C9 = 5950 x .392 = 2330 
U 
a = B20: oie T, = 2200 - 965° 760 _ 1345°R 

O 6160 ar Sol “Feary: 

Va _ 5950 
= ‘ I 64.0 M = 2 = —_—_—_—_ = 2.2 

An 72.0 N1345 = 2 D a 2640 y) 


3810 
SSS egal 
Mp 2610 a 





Site 


VR Wo _.91 x 3810 


Estimate M, = ~ hmmm ————— elon, clin Wye ely (Shaeyes 3} 5 
oe 8) owt, 720,05.) pe ee 


(c) Determination of rotor velocity coefficient, Up 


ay = .48 sin® 60° = ee) (Yip = 8 x .03 sin© 60° = alts 


ase 


(Y ~ .12, and estimated Oe 


PR TR esas 


From Table I-d, for Many a ae. Vp = neil ee 
mam estimated Vp = .881 is obtained for comparison from the subsonic 


reference curve of Figure 6. 


(d) Completion of velocity triangle and performance calculations 





ee ort X 5010 = 5520 ft/sec ee = VeVax,2 = .O71 x 1900 = 1655 
V, = V1655" + (3320 cos 60° - 2330)° a = cos) (2022) - 18.2° 
1740 
= 1740 
mean radius, R= 60 x 2330 x 12 = lLyeetel 
ex x 15000 a 
: O 
ei Pee WwW = 2330 x 3810 x 1.871 sin 60 ee ae 
245.36, Le 
ae 275 _ 756 GO _ f0T x 28800 _ 35.4 1b/sec 
Msg) © S10) UE 5D saan 


The calculated efficiency appears reasonable but perhaps 


eee cly optimistic, compared with Figure 1. 


(e) Temperatures and total head efficiency 


2 2 
. Wo 2) _ 3810 Z O 
Ty = To =P DeJICp (1 - Vp ) = 1345 a 4.13 x 10" x enn = Wee R 





zs 2 
- Vi - 1740 — 1502°R 
42a. UO + ge 





aS 


W = CpAT! = .827 (2200 - 1502) = 577_Btu/1b. 
Wh 3320 
=n Sr ee 
4” 2.0m, 2.0mm 
Pe 20 = agi Sa a 
72.0 V0, 72.0N 1429 se 
Pts PB, (1 + Zo 2 50 (1 + .125 x .642)9*° = =8.5 psie 
W 575 
= a = —] = = Ro. 
"r = An. * W827 x e200 [2 - (ey ey 
Is 450 


Nozzle Flow Analysis and Profile Layout 


A detailed analysis of the stator losses and required geometric 


layout is necessary. Pertinent assumptions are: 


(a) 
(b) 


flow within the nozzles (stations O-l1) is considered isentropic, 
all stator losses occur in and downstream of the exit plane, 
between stations 1] and 2, and a uniform flow is established 
beyond station 2. The losses between 1 and 2 are represented 
by the coefficient y, = O65. 

Scibinuicy 1S satisfied beuween lL anc 2 

tangential momentum is conserved between 1 and 2. 

the flow is adiabatic, Te = TU: 


The continuity equation is developed analogous to that 


expressed by Equation (II-6-7) except using total temperatures and 


pressures. The resulting relation is 


M evh. 
cos Oy = &). + Cp ea (b+ 125 My") ? 


Cos 5, 
Mi ieee My@ 





aie 


Prom previous calculations, On = false M, = 2.2). The pressure 
Seooitticient , 
ce ee Oi Ole 25) x Dp 
p= =f ae eee = CTT 
Ee 450 


feamec. there is 


2\4.5 
a (1 + .125 M,“) 
cos @, = Ge): + .0613 a (2) 





Tangential momentum conservation gives Vy, sin On = Vo sin An » 









M, VYgRT, sin a, = M, V7gR am sin Q, = Mp J72R te sind, 
eee 1+2/ 
2 2 
or 
Yl 2 JUs .125Me 
sin Q = 2 eel ras a Sin Q, = ee -<1-> Vit -125MT sin (plow 
1 14 y=-1 M2 Nay 25) 3 2.25" My 
5 2 
V1 + .125 Me 
sin Q) = 1.6/0). == (2) 


My 


For a nozzle with no trailing edge thickness, simultaneous 
solution of (1) and (2) gives M, = 2.150, @ = 77.2° . 
For structural integrity a nozzle having 10% "blockage," or 


t/S cos A, =.10, is selected. Equation (1) then becomes 


(1 4.225 M2)" a 


cos Oy = .0682 





M) 


Solution of (2) and (3) gives OQ, = 7505" M, = 2 7o hor 





the required nozzle exit conditions consistent with the desired (velocity 
triangle of Figure 12) discharge flow. The blade angle turns out to be 


about the largest practically acceptable value. If desired, the magnitude 





at 


of & could be decreased by providing a slight annular expansion 
Bepween Stator and rotor. 

The type of nozzle profile is probabil, notycrauleaine 
"sharp-corner" nozzle of the type discussed in Reference 16 was chosen 
somewhat arbitrarily. The sharp-corner nozzle, in general, provides 
for the shortest possible nozzle length capable of a uniform discharge 
flow. Its properties become much more favorable in cases having larger 
nozzle pressure ratios and Lower discharge angles than the present 
example. 

Rather than using a full characteristics network, the approx- 
imate nozzle contour required was obtained by calculating the required 
bmiroat angle, over-all area ratio, and the approximate orientation of 
ioe limiting characteristics. The angle of turn at the throat is 
simply = where v, = v3 (My) the Prandtl-Meyer property angle corre-— 


Bpending to the exit. 


i 3 tan7+ | i (M, = - 1) 2 Hansen 


D 178 ==1T- 9, the required angle. The requared area 


Vv 


ek 


facto is obtained from one-dimensional flow theory. The relation is, 


(e.g. Reference 17, problem 4.20), 


Pi Val ey 
as (—-) / 


Py 2 yar fa(y-1) 
YT 2 1 \r-4/7 Ca 
t | ae iC : =| y+1 


Beem, = 2178, 


IL 


ae 
yee (le - 12> x 2.1782)? = 10.22 





nee 





and 
A - FO 45 
_ 7 len 22)) (8) _ 2.90 
t V8 (10.2240-1 9 
tne schematic nozzle layout as showneinePicire ss 
The limiting characteristic of the throat fan, ab, stands 
. 7 O ee dl: ae O Pe iy | ie 
ou approximately p +u=- 1l/.9 + Sin [ 7 alg e =) eo cee Gman= 


me.6. The limiting reflected characteristic, bce, stands at approximately 


es) = 27.4°. Thus the minimum nozzle length is 


ie Sin 


We 2-20 Wt 


pe Cee = (OB ee = 7.0 
tam 19,6° iY Heya eeg ue ( if D) oh OD oh 


fewever, due to the required truncation of the nezzle at =" Q, = aig race 
mioeaGimension be must be no less than ae0 Wt = 8.38 W,+ Hence, 


ten 172 
minimum throat to discharge dimension, parallel to nozzle axis, is 


6b +bg = 2.0 + 0.4 wm. = 11.2 w,. The resulting nozzle profile layout 
is shown in Figure l4a. 
Since the edge thickness, t, is 1/10 the total width, 


t = 1/9 x 2.20 wy = .2hh w,. Total width is w = 2.444 w,. ‘The spacing 


i OM we 


between blades is s = cos Gy. = Gos 75. 3° = 9.65 We The number of 


blades Z, = 2@R = en xl -9 _ 11.59 . For a selected throat dimen- 
Ss 9.05 Wt Wy 

sem Ww = .290", Z, = 40 blades. The nozzle cascade is illustrated 

in Figure lb. 


The height of the nozzles (and annulus), may be more pre- 


Giocly calculated, 


oy 
hy - WR To = 25 ae 12070 x L355 googie fo 
ee any 2 emk 40 X 1900 X 2x x ILE fe. ee 


The total throat area is Zghw, = 40 x .961 x .290 = 11.14 in® . 
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With the nozzle dimensions fully determined, the calculated 
stator coefficient of III(a) should be checked. For a thickness ratio 
t/c meri and Solidity c/s = 1.44, Figure 2a and Equation II-2-1 gives 


(Yp) 5 = = .05 = .053. The secondary losses are (Y 


secs 
ne O 

-15~= sin® (0-3 = .272, where the blade angle rather than the gas 

.961 cos 62.3° 

angle is employed to provide the most conservative estimate. Total 

losses igaee ta re = .33 exceed the estimated .28 of previous calculations. 


Hence, from Table I-d, the estimated ¥, = .965 may be about .004 too 


optimistic, but recalculation of IV-4 seems unnecessary. 


5. Rotor Flow Analysis and Profile Layout 

A leading edge thickness ratio, (yy) Se must be selected first 
so that incidence may be determined. For a chord of about 1", an 
expected solidity c/s = 2.5, and a desired edge thickness of about 
Meee there is t/s = t/c - c/s s .020 x 2.5 = .050- 

The incidence, Bo - Bz may be determined from Equation 


cos 
= + —- . From considerations discussed in 


(II-6-10), cos Bz = 
II-6, it is assumed that about one-quarter of the rotor losses occur 
around the leading edge - between stations e and 3. 
For (Yrn)p = Lie = 66 _ 165 there results f 
ile LE/R ~ 7 ib oe one . rom 
Mele I-c, (Cp)rp = Cp(M = 2.44, Y = .165) = .90 


O 
cos 60°) = 52.8° and the incidence is ae ce 


Hence Bz = cea) (Gone = 
The assumed pressure loss has increased the incidence by about Sch 
from the isentropic case. It is of interest that, were the blade 
engle, Pa, equal to 60°, the relative incident gas angle would be 


about 66.1° and the nozzle pressure ratio would be greatly reduced 


from the desired 15 to 1 value. 
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The rotor profile was laid out following the scheme of 
Reference 14 as summarized in Section II-5 of the present report. 
As discussed in Section II-6, the rotor entrance Mach number, Mr? is 
assumed equal to Min = 1.44. It may exceed this value slightly, as 
suggested by the discussion of II-6. 


By Equation II-5-1, 


nae 
(eran 2 208 Hs = cos (sin —— ) 

out 3 

rs - ay 
—— > cos (sin7t ——) - : (20 
ee ren ma ( Te ‘ 
Ey _ — ON ee eee ee 8g 
ry’ bran Vacs ae lL + .720 


@ 


Theoretically a channel width-to-radius ratio of .426 would 
Just insure tangency of Mach waves with the convex surface, while allow- 
ing an excessively large Prandtl-Meyer acceleration on this surface. 
Hence, the authors recommend at least a 20 to 30% reduction in Galo 
Accordingly, a design value, bran = .22 was selected. For an 
waeecary construction scale, let (ryjiren = 1-0s Wiman = Ween = ees 
(Po)tran = 1-11, (ri )t+an = -89- Since the curvature is to be doubled 
in the central section, (Ty) een = 4 50), nese eae pene le CAP. = eck 

The profile is constructed, to arbitrary scale of 1.10" 
channel width, in Figure 15. The transition curvature is 15° and the 
leading edge consists of a WO wedge. This is slightly less than the 


maximum angle consistent with an attached shock. The edges could prob- 


abeebe rounded at a slight sacrifice in performance. From construction, 
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the solidity, c/s = 2.88 and the thickness ratio t/e = .160. For a 
selected chord of 1", thickness at the shoulder chamber is approximately 
mezO" and the scale of Figure 15 is 5.74/1. 


The number of blades required is calculated from 


i Get ERR 2 Bae Bn eee ede 
S Ss. T SSS = 
Gee ewes} 


The rotor losses may be recalculated more precisely for com- 


parison with (Yn)p = .66 obtained in IV-3. For t/e = .16 and s/e = .35, 
Figure 2b and Equation (II-2-1) gives (Yp) . 26 x23 = .104. The 


R 
revised best estimate of rotor loss is (Yp)p = -64. From Table I-d, 


Vp = VWp(Y = .64, M = 1.22) .668 , 


compared to the previous estimate of .871. From Equation (III-2-3), for 
the final estimates of ¥,, Vp, the final efficiency estimate is y = .748. 


From IV-4, hp = hz = 7 Ooh a amel 


a = A ae sae (961 Lie? 2 eae 
To VR 1345 .868 
From IV-2, 6 = 1910 Rh, = 1910 x 17.8 x 1.172 = 39800 psi 





B = ony, > cos B = 39800 x .748 x .868 x .64 cos 60° = 8270 





2 
Total weight, m==3220 »° = 3220 (.392)" = 495 lb. 


6. Discussion of Resulting Design 
The major numerical results are compiled for comparison with 


estimates of IV-2. 





=e 


n = .7h8 B= 86270 ac 
ae 35.4 lb/sec. B = 60° 

cg = 39800 psi 7 = 1.44 

Ree = ky Or My, = 1.20 

¥, = -961 9 = .392 

Vp = .868 6 = .6h 


The calculated efficiency of .748 agrees well with that 
indicated in Figure 1 for Wc = .478. Hence, a confident prediction 
of efficiencies well above .7 offers encouragement for improvement 
over the previously demonstrated performances of supersonic turbines. 

ihe Design Point Selection Vechnique, leadimsere siiel yee 
proves to be of great assistance in presenting a thumbnail sketch of 
meewlable design choices. The slight variation from the predicted 
values of IV-2 was obviously due to the reductions in both Vp and Ve 
from the approximations assumed. Of particular interest is the observa- 
tion that the "invariant" B changed only .24% from the value based on 
the original arbitrary design point. 

Perhaps the most questionable area of the problem was that 
Seeestimating rotor incidence. The theoretical basis is certainly 


open to speculation, as discussed in II-6. 





V. CONCLUSIONS AND DISCUSSION 


This investigation indicates that well designed supersonic- 
rotor single-stage impulse turbines are capable of rather impressive 
performance. These machines appear to be an attractive compromise 
Beoween simplicity and high efficiency. Utilizing recent test results, 
an effort is made to better understand the losses and flow peculiarities 
ewoupersonic stages. Many questions are encountered which cannet be 
@ecisively resolved without more tests of the quality previously reported 
by the Aeronautical Research Council of Great Britain. 

Analysis of supersonic turbine tests indicate that no appre- 
Muele Losses need be attributed solely to Mach number eifreets.  pro— 
famed cCerbain design criteria are followed. The major rotor blade 
Seater Objective is the prevention of shock-induced flow separation. 
fretiles have been recently developed by the Aeronautical Research 
Seunet] which virtually eliminate separation. In supersonic flow these 
blades have much lower profile losses than the conventional constant- 
curvature blades commonly used in subsonic rotors. Carefully designed 
profiles yield considerably greater improvements in total efficiency 
mem static efficiency. Properly chosen rotor blade flare, accurately 
mabeeting actual losses, has a considerable effect on static efriciency. 

Since Mach number level need have no direct effect on blade 
losses, an extension of the Ainley-Mathieson loss prediction procedures 
to the supersonic case is proposed. Some straightforward modifications 
to the basic procedure are adopted. The extension to supersonic flow 
appears tentatively justified by the meager amount of available test 


results. 


ee 
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The existence of a unique rotor incidence, in supersonic 
flow through blades of finite edge thickness, profoundly influences 
pressure and velocity distributions throughout the stage. The 
incidence results from continuity requirements which demand a turning, 
mac probably an increase in velocity, at the rotor entrance. ‘The 
explanation of the phenomenon presented in Reference 14 and based on 
Prandtl-Meyer corner flow is shown to be unsatisfactory. An alterative 
explanation is developed based on the assumption of negligible velocity 
increase. More experimental work is required to yield a clear under- 
standing of this important phenomenon. 

A technique is developed to assist the designer in selecting 
the "optimum" design point for a single-stage impulse turbine having 
specified RPM, power, and blade root stress. This method is independ- 
ent of Mach number level and the method of estimating blade losses. 

It greatly simplifies the decision by promptly resolving design point 
Pebection criteria into a readily digested form. 

The preliminary design of a large rocket-turbine is developed 
Bena numerical illustration of the foregoing procedures. Blade speeds, 
but not necessarily stresses, are typically very high in order to 
reduce leaving losses. The results indicate static efficiencies 
exceeding .7 may reasonably be anticipated from well designed turbines 


Sieonis type. 
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APPENDIX A 


NOMENCLATURE 


Gige a. 

function defined by Equation III-4-11 
chord 

Constant pressure Sspeciti1e cae 


pressure COEiTiclient; ratio Of eectuel to isentropic tovak 
pressure 


isentropic jet velocity corresponding to stage pressure 
ratio 


function defined by Equation III-4-12 
hOonee 

function defined by Equation III-4-lha 
function defined by Fquation III-4-14b 
gravitational constant 

enthalpy 

Dikade or annulus heieat 

rotor tip clearance ratio 

oi ten Mc Lene es a nailie 

mechanical equivalent of heat, 778 Btu/1b. 
Mach number 

mass 

RPM 

pressure 

mean blade radius 


gas constant 


Lies 





BGs 


S - blade spacing 
Ay - temperature 
16 - blade edge thickness 
U - mean blade speed 
V - absolute gas velocity 
W = gas velocity measured relative to moving rotor 
W - mass flow rate 
al - Specific work 
W - channel width 
loss parameter. 
Z - number of blades 


GREEK LETTERS 


ee - absolute gas angle 
B - relative gas angle 
Y - ravio Of Specific heats 
= Wo/Vp 
Se Hig - Myo 
TW 5 
9 = U/Ve 
nN - Seconcaary Loss #2cCelergurenesparameuet =o! Higure 4. 
p - wall angle 
V - velocity coetficient; ratio of Setual to isentropic 
velocity 
i = efficiency (static) 
0 - COLor blade Ccenuriiusal root Stress 


O - density 


— 





maximum blade thickness 


Prandtl —Meyer angle 


sin-t @: the Mach angle 


enpular velociuy 


SUPERSCRIPTS AND SUBSCRIPTS 


total gas properties referenced to absolute gas velocity 


total relative gas properties) renerieed. tem row emecic 
enlore al wp's 


sesimene sakedlicne 

Staver Gxt) plane 

downstream of stator, after mixing 
HUStealver Toler Teadangeedeas 
downstream of rotor, after mixing 
Protrile-4 Of Section Wiles 

axial component 

profile B of Section II-5 

cross section 

Cenural or Mid-—chord Toler cnanne 
inner or convex rotor profile surface 
LSehicrop Te 

rotor tap loss 

hebor leading edge 

MeamW WOver jprotl leveunvature 
vector mean angle 


outer or concave rotor profile surface 





profile loss 

rotor 

stator 

secondary loss 

total loss, or total head efficiency 


nozzle throat 


- rotor entrance transition channel 


COPreSpondi nes toyrelative. Vell@ein, 
UPS tral, Or s1mleG, peo tren 


downstream, or exit, position 





APPENDIX B 


DEVELOPMENT OF EQUATION II-3-1 


By reference to a T-S diagram one may easily show (e.g. 
Appendix B of Ref. 4) that 


BIL ey" _ PL'/Po' 
ea 
PS / Po Cua) oo) <a ameen (ce) ean ci 


(Cp)p y/ y= 


t + { t 
ty, ae may now be expressed in terms of qn and p), /P., 


With Ty, ' 19? the temperature corresponding to an isentropic pressure 
3 


deleranag -1 
decrease to p)', there is Ty, Is = 1," (P,*/P,1) 7 7 and 
(y-1)/7 
Vie ' a ' _ t t 

Lo = ae age = ue [1 (P), i ) ] 
Since 

; - ee 2 ie _ te _ tt 

ie i Pee eee! tT! a (P '/P ry Y-L)/ 7) 

O 4,19 O Bene lea. 

there is 


NY 
10) EGR St ae tee en 5) 


and finally, 


(Cp)R = Seno TS 
(Co), [ (1-1) a3 nn (Py, '/Pot) 7-2/7} y/y-1 
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Figure 1. Comparison of Theoretical Efficiency Estimates for 
(Subsonic) Single Stage Impulse Turbines with Results 
of Supersonic Turbine Tests. 
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Figure 2. 
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Figure 6. Estimate of Rotor Velocity Coefficient - 
Typical Impulse Blade in Supersonic Flow. 
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Figure 7. Cascade Loss Results of Reference 13. 
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Figure 8. Comparison of Rotor Incidence From Experiment 
and Two Theories - For Turbine of Reference e. 
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Figure 9. Velocity Triangle Nomenclature. 
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Figure 11. 
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Constant Stress Efficiency Curves for Turbine 
Design Problem. 
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Figure 12. Design Velocity Triangle - 
Mean Radius Conditions. 
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“igure 13. Sharp-Corner Nozzle Geometry Considerations. 
(Not to Scale) 
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Figure lia. Design Nozzle Profile. 





Wy =.290° 
tT/e =.2\ 
c/s = 1.44 


SCALE: 1/1! 


QO 
So 230" 


ee 


Figure lib. Sharp-Corner Nozzle Cascade 


















